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Abstract 

Toponium production at future hadron colliders is investigated. Per- 
turbative QCD corrections to the production cross section for gluon 
fusion are calculated as well as the contributions from gluon-quark and 
quark- ant iquark collisions to the total cross section. The dependence 
on the renormalization and factorization scales and on the choice of 
the parton distribution functions is explored. QCD corrections to the 
branching ratio of rj t into 77 are included and the two-loop QCD po- 
tential is used to predict the wave function at the origin. The branch- 
ing ratio of rjt into jZ, ZZ, HZ and WW is compared with the 77 
channel. 
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1 Introduction 



The production of heavy quarks at hadron colliders has received a lot of attention 
from the experimental as well as from the theoretical side. Higher order QCD cal- 
culations for open top quark production |], Q provide the theoretical basis for 
the current experimental limits 0] on the mass of the top quark. Given the cur- 
rent range for the top mass of 90 to 180 GeV deduced from electroweak radiative 
corrections ongoing experimental studies at the TEVATRON should be able 
to discover the top quark in the near future. 

Planned experiments at LHC or SSC will lead to a huge number of ti events 
and may pin down the top quark mass with a remaining uncertainty of perhaps 
5 GeV pi . Aside from the intrinsic interest in the precise top mass determination, 
the aim to fix the input in one-loop electroweak corrections |F] constitutes one of 
the important motivations for the drive for future increased precision in the top 
mass. 

In H it has been proposed to decrease the experimental error significantly 
through the study of rjt in its 77 decay channel. Making use of the large luminosity 
of future hadron colliders, one should be able to overcome the tiny cross section 
multiplied by the small branching ratio into two photons. The excellent gamma 
energy resolution, originally developed for the search for a light Higgs boson, 
would in principle allow the measurement of the mass of the bound state at 
about 100 MeV accuracy. This would result in a determination of the top quark 
mass limited only by theoretical uncertainties. It was shown that the toponium 
signal could be extracted from the 77 background with a signal to noise ratio 
decreasing from 2/3 down to 1/8 for m t between 90 and 120 GeV. 

The estimates of || for the production and decay rate were based on the Born 
approximation and on fairly crude assumptions on the bound state wave function 
which was derived from a Coulomb potential, neglecting QCD corrections. The 
signal to noise ratio and the statistical significance of the signal as estimated in 
|| is just sufficient for the discovery. Therefore a more precise evaluation of the 
production rate is mandatory. 

In lowest order the reaction is induced by gluon-gluon fusion (fig. la). QCD 
corrections will be calculated in this work. They affect and indeed increase the 
rate in this channel by about 50%. Quark-gluon scattering and quark- antiquark 
annihilation into i] t l enter at order a^. As we shall see the qg process increases 
the production cross section by about 8% whereas the effect of qq annihilation is 
negligible. 

The production cross section for r] t with nonvanishing transverse momentum 
is infrared finite in the order considered here. The total cross section, however, 
exhibits singularities that have to be absorbed in appropriately chosen and defined 
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parton distribution functions. The dependence on the factorization as well as 
renormalization scales will be studied as well as the dependence on the choice on 
different parametrizations of parton distribution functions. A full NLO calculation 
requires the corresponding corrections for the physically interesting branching 
ratio. These are available and will be implemented. 

The bound state production cross section and the branching ratio of rj t into 
77 depend (markedly!) on -R(O) 2 , the square of the bound state wave function 
at the origin and hence on the potential. The range of -R(O) compatible with the 
current knowledge of the perturbative two-loop QCD potential and the present 
experimental results for Aq C d will be explored. 

As we shall see in the following, fairly optimistic assumptions are necessary 
to open even a narrow window for toponium discovery in the mass range of 
m t = 90 — 110 GeV. This is the consequence of the dominance single quark 
decays (SQD) and the small branching ratio into two photons. The situation 
improves dramatically if new hypothetical quarks are considered with suppressed 
single quark decay. Examples are b' or isosinglet quarks with small mixing with 
ordinary d, s or b quarks. These would dominantly decay into gauge and Higgs 
bosons, leading to spectacular signatures. 

The outlay of this paper is as follows: In section 2 the QCD corrections to gluon 
gluon fusion into r\ t and the cross section for the qq and qG initiated processes will 
be derived. In section 3 the formulae for the decay rates of r\ t will be listed and the 
dependence of the wave function on the choice of the potential studied. Section 4 
contains the numerical evaluation of the production cross section and the study 
of scale dependencies. Section 5 contains our conclusions. A brief account of this 
work has been presented in M. 

2 QCD corrections to hadronic rj t production 

For the complete evaluation of NLO corrections for r\ t production in gluon gluon 
fusion virtual corrections (fig. lb) are required as well as corrections from real 
gluon radiation (fig. lc,d). In addition one needs the cross section for quark- 
antiquark annihilation into r/j+gluon (fig. 2a) and the (anti-)quark-gluon (fig. 2b) 
initiated reaction into r^+quark. Technically it is quite convenient to employ 
dimensional regularization for both ultraviolet and infrared singularities. The 
calculation of amplitudes for bound state production (or decay) can be performed 
in two different ways. One may either evaluate directly the amplitude for the 
production of a (QQ) bound state with the desired spin and orbital angular 
momentum configuration and arrive at an amplitude proportional to the wave 
function at the origin or its derivative ||T0| , |TT| ; or alternatively one may evaluate 
the production rate for open QQ and subsequently identify the rate at threshold 
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with the rate for bound state production, substituting 

dPS n (ki + fa; fa,..., k n , k Q , kg) — > dPS n ^i(k x + fa; fa,..., k n , q) (1) 

where M = 2m t + E Bind and kg = kg = q/2. For the one particle final state 
relevant for resonant rjt production in gluon gluon fusion this implies (in 4 — 2e 
dimensions) : 

dPS 2 (fa + fa; k Q , kg) ^2tt5( S - M 2 ) (2) 
and for the r] t + parton configuration that will be of main concern in the following 

dPS 3 (fa + fa;fa,k Q ,k ) => J ;dt (3) 

where the Mandelstam variables are defined by 

s = (fa + fa) 2 t = {fa- q) 2 u= (fa-q) 2 = M 2 -t-s (4) 

In the latter approach it is important to make sure that — in the limit of vanishing 
relative velocity — the cross section receives solely contributions from the QQ spin 
parity configuration corresponding to the bound state. It is therefore mandatory 
to restrict the evaluation of the amplitude to QQ color singlet configurations and 
in the reaction under consideration to those three gluon states that are totally 
antisymmetric with respect to color as well as to their Lorentz structure. 

The former approach leads typically to more compact expressions^. However, 
for states with J PC = h the formula for the bound state amplitude involves 
the 7 5 matrix whose formulation in the context of dimensional regularization 
introduces complications. This purely formal difficulty is absent in the second 
approach, which for this reason has been adopted in this calculation. 

2.1 The hadronic cross section 

In this section we define the general structure of the cross sections for rj t produc- 
tion in hadronic collisions 

h 1 (P 1 )+h 2 (P 2 )^ Vt + X ^-f + j +X (5) 

within the framework of perturbative QCD. Here fa, h 2 are unpolarized hadrons 
with momenta Pi. The hadronic cross section in NLO is thus given by 

a H (S) = J dx l dx 2 ^K^Ql)^^2,Q 2 F )a a \s = x 1 x 2 S,a s (ii 2 ),ii 2 ,Q 2 F ) (6) 



2 See for instance |fl3| for particularly impressive examples. 
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where one sums over a, b = q,q,g. fa(x,Q 2 F ) is the probability density to find 
parton a with fraction x in hadron h if it is probed at scale Q 2 F and a ab denotes 
the parton cross section for the process 

a(ki = x x P x ) + b(k 2 = x 2 P 2 ) ^rj t + X (7) 

from which collinear initial state singularities have been factorized out at a scale 
Q F and implicitly included in the scale-dependent parton densities f£ (x, Qp). We 
next specify the possible partonic subprocesses that contribute to r\ t production 
in LO and NLO. 



2.2 Gluon fusion 

In Born approximation the cross section for gg —>■ rjt (fig. la) is given by (Nc = 3) 
with 

AT 2 - 1 

\MBam\ 2 = "g^-fli 16 i 1 ~ ^X 1 ~ e ) ( 9 ) 

The QQ state has been projected onto the color singlet configuration. The division 
by 64 • 4 • (1 — e) 2 is implied by the average over 8 2 gluon colors and (n — 2) 2 gluon 
helicities in n dimensions. g s will be replaced by a s through g 2 = 47ra s /z 2e . The 
Born cross section can therefore be cast into the form 

where z = M 2 /s. Virtual corrections to the rate for QQ production close to 
threshold (fig. lb) can be combined to yieldQ []i~2] 

/47m 2 V 

virtual Born Born I -TFT r (! + e ) 

a s 7r 2 , 1 , T 1 ,1 ,\ 

[Cp— + b N c — -bo— + A) (11) 

n \ v e uv ej R e IR ) 



where 



A = -T F Un2 + C F (^-^+N c (^K 2 + 
b = Nc^r-n f T F ^ (12) 



3 Here C F = 4/3; N c = 3; and T F =l/2. 
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1/ejjv and 1/ejn represent poles arising from ultraviolet and infrared divergent 
integrals respectively. The mass singularity 21n(rrtj/M) in |L2] has been replaced 
by 1/ejR. Since we are considering bound state production the 1/v term (v de- 
notes the relative velocity of the two heavy quarks) in eq. (|Tl"| ) is absorbed in 
the instantaneous potential and thus effectively dropped. The QCD coupling is 
renormalized in the MS scheme and is given by 



a s _ a MS (^2\ 



l-^b (—+lu47T- lE 

vr \e uv 



(13) 



7T 7T 

Combining Born result and virtual corrections one thus obtains for resonant pro- 
duction 



" gg / \ 

° Born+virtual ~ a 0{ e ) 



1 + *(*£) ni + t 
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-N c — - b — + A 



1 h 1 

2 b 0~ 

e IR e IR 

where the normalization factor 

, l7r 2 i? 2 (0) 2 , 2 , 4e (l-2e) 

has been split off for convenience. 

In a second step the cross section for 



S(l-z) (14) 



gg^gQQ (16) 

has to be calculated (fig. lc,d). Only the kinematical situation will be treated 
where the relative momentum of Q and Q vanishes. In an intermediate step the 
matrix element for 

999 ^QQ (17) 

will be considered (fig. le,f) and reaction flTBj ) can be subsequently obtained 
through crossing. The QQ state will be projected onto a color singlet configura- 
tion through 5ij/ y/3. The overall wave function of the gluons is of course totally 
symmetric. Requiring a totally antisymmetric configuration in color space im- 
plies, in turn, a totally antisymmetric wave function in Minkowski space. This 
latter condition secures J PC = CT + for the QQ bound state. (Conversely, requir- 
ing symmetric wave functions in both color and Minkowski space would lead to 
J PC = I™.) 

The amplitude arising from the diagram fig. le is given by: 
A{k 1 ,k 2 ,k 3 ;e 1 ,e 2 ,e 3 ) = 

u m + k2 + y/2) H9 M f_A^+M^ ( h ) 

[{—hi + k 2 + «3j /4 - m 2 \ 
A- M+ /k 3 )/2 + m] . , 1a x 
" [K-fe + W-m^] { ~ 19 43> {{h + h + h)/2) (18) 
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to be multiplied the color factor 

/ A ai A a2 A a3 \ 5 - 11 

(tttI^^Iw-^' (19) 

As discussed above, only the / term will be retained. The amplitude is obtained 
from the six permutations of the diagram without triple boson coupling 

i 1 

«^ nr ^ f a\a^a-i ^ ] Qjm-^^ij kjf km-, ^ii ^ji ^m) (20) 



The second type of amplitues, represented by fig. If, involve the triple gluon 
vertex. Their sum reads as follows 



B(h, k 2 ,h;ex, e 2 ,e 3 ) = 
u ((A* + A; 2 + fc 3 )/2) 



z [(- fo+ fc 3 )/2 + m] . 



(fci + fc 2 - «3j /4 - m 2 



i; ((A* + k 2 + A; 3 )/2) (21) 
x ,, 7' , 2 (-p) {(/ci - A; 2 ) a ei e 2 + (A* + 2/c 2 )e 1 e 2 t + (-2^ - fc 2 )e?e 2 } 

(fci + k 2 y 1 - J 



with a color factor 

/A*A-\ ^ _J_1 

^22 y ^3 26 v/32 Jaia2a3 1 J 

Three cyclic permutations are combined in this class of diagrams 

i 1 2 ~ 

^ ~~ To 4 f a i a 2 a 3 ^ / ~r^{kii kj, k m , 6j, 6j, 6 m ) (23) 

V i,j,m(cyclic) 

The ghost amplitude C is easily obtained from eq. fl2~T| ) through the replace- 
ment of the curly bracket {. . .} — > k 2 , where k\, k 2 and fc 3 denote the momenta 
of the incoming ghost, antighost and gluon respectively. 

The squared matrix element is finally obtained in n = 4 — 2e dimensions 

EM* = El*4 + 3| 2 -2EKf = E iC 2 ^ 6 128(1-6) (24) 



st + tu + us 
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(s-M 2 )(t-M 2 )(w-M 2 ) 



48 

01,02,03 ^ U 

M 8 + g 4 + t 4 + w 4 t 
stu 



[1 - 2e) +4eM 



where 

s = (/c 1 + /c 2 ) 2 t = (^i + fcg) 2 u= (A;i + & 3 ) 2 ; M = 2m + E Bind (25) 
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After crossing (fc 3 — > —k 3 ) this result is applicable to the reaction of interest 
(Hp. In the QQ threshold region (with the color restriction discussed above) one 
obtains 

da(gg -> gQQ) = - - _ £ \M\ 2 R dPS 3 (26) 



■■=9 6 \M\ 2 R 

where gluon helicities and colors were averaged as usual. Replacing the three- by 
the two-particle phase space as indicated in eqs. flip and (|p one obtains 

d°{gg -> grjt) = 00(e) — jj — (27) 
/ st + tu + us \ 2 M 8 + s 4 + t 4 + u 4 / 4ns fi 2 \ e dt 



V(s-M 2 )(t-M 2 )(M-M 2 )/ stu V ut J r (!- e ) 

The term proportional 4eM 2 in eq. (|24]) without singular denominator will not 
contribute in the limit e — > and has therefore been dropped. In the limit e — > 
the result coincides with eq. (8.45) of [|13|]. For the subsequent discussion it will be 
convenient to follow closely the treatment of NLO corrections for the Drell-Yan 
process presented in the classical paper by G.Altarelli et al. [fL4| . The invariants 
s, t and u are expressed in terms of M 2 , y and z by 

M 2 M 2 M 2 
s = —; t = I-* U = (l- z ) y (28) 

z z z 

and one arrives at Gluon 

°{gg^g{QQ)) = * (e)-r(i + e J lrr/r ^ 

7T 



6 r 1 



ri+eri 



m 2 J 

[ L f 9g zf(l-z) 1 - 2e y- e O.-y)-<dy (29) 



The normalization factors have been chosen to allow for easy combination with 
the virtual corrections and f" is given by 



z ( (1 - z)(l - y)y 



4 \(l - Z )[l - y(l - Z )][(l - Z )y + Z ] 



X 



z* + 1 + (1 - ^) 4 ((1 - yf + y 4 )} ( - + (30) 



As a consequence of the symmetry of f" with respect to y *->■ (1 — y) it suffices 
to evaluate the 1/y term. The y integration can be performed with the help of 
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the distribution identities 



x 



-1-e 



z e (l - z)- l - 2f - 
y- 1 - e (l-y)- e 



--5(x) + (- 



\nx 



x 



In z 

+ e- 2e 



1 - z 



'ln(l-z) 



c ln(l - y) 



(31) 



and 



1 



,7T 



1 - t z — + . . . 



r(i-e)r(i + e) * " 6 1 "' (32) 

Combining the Born result with real radiation and virtual corrections one obtains 



a(gg^rit + X) = a (e) {5(1 - z) + 



a s ( An fi 2 



71 



M 2 



r(i + e) 



b In j£ - iV C y + ^ 5(1 - a) - ~P g9 (z) + iV c F(^ 

where &o an d A are given in eqs. (O) and 

ll^ 5 + lLz 4 + 13z 3 + 19^ 2 + 6,2-12 



(33) 



F(z) 



0(1-*) 



+ 4 (- + *(! -2 ln(l-z) + 4 



6^(l + 2) 2 

'ln(l - z)\ 



1-2 



(34) 



+ 



'2(z 3 - 2z 2 - 3z - 2)(z 3 -z + 2)z ln(z) 



(1 + ^) 3 (1 -Z) 



1 - Z 



and 



P n = 2N c 



+ 



1-z 



+ z(l-z)-2) + b S{l-z) 



(35) 



The 1/e term can be absorbed in the gluon densities. If these are defined at a 
factorization scale Q 2 F , one effectively subtracts from the above result the term 



^o(e)^ 



a s ( 4:71 yU 2 



2tt V Ql 



r(l + e)[-~P gg (z)+C gg (z) 



(36) 



In the MS scheme C gg = 0, in the DIS scheme 



C 99 (z) 



-n f 



(z 2 + (l-z) 2 ) l n ( — )+8z(l-*)-l 



(37) 
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The partonic cross section is then given by 



a 



uu 



1 7T 2 R(oy 



«^ 2 ) *(i - *) 



+ 



s 3 M 3 ~ MS { 



~ ^P 99 (z)-C 99 (z) + N c F(z) 



TV 



C F — 



3^/111(2) 



(38) 



2.3 Gluon quark scattering 

Gluon quark scattering g(ki)+q(k 2 ) — * q(k 3 ) + r] t (P) proceeds evidently in leading 
order through the diagrams in fig. 2b. Collinear singularities have to be ab- 
sorbed in the structure functions. The (dimensionally regularized) squared matrix 
element for the production of QQ at threshold is given by 



\M 9q \ 2 (s,t,u) = 32 



s 2 + u 2 1 . . e 
1 -2e) + - 



(39) 



and has to be supplemented by the color factor 



E 



2 LI 2, 



A7 



1 

3 



(40) 



and the coupling constant g®. After spin and color averaging the cross section can 
then be cast into a form identical to eq. ( p9f) with /" replaced by 



/ 
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1 



l+y\l-z) 



[l-y(l-z)y l-2e 



9(1-*) (l-y) 

The j/ integration is performed as before and only a single pole remains 

,2\ e 



(41) 
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99 



a ( e )^r(l + e) 



47r/r 

1^2" 



+ In 



P 99 (z) + C f 2(1 ^ ^ (1 - \nz) + C F z 



(42) 



with 



c, 



l + (l-z) 



(43) 



The 1/e term can again be absorbed by the (quark) structure function. This 
amounts effecively to the subtraction of 



-P gq {z) + C gq {z) 



(44) 
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from the above result. The function C gq is again specific to the subtraction scheme: 
C gq vanishes is the MS scheme and in the DIS scheme it is given by 



C gq (z) 



-C m {z) 



(45) 



where 



C qq (z) 



C f 



:i+z 



2 . An(l-z) 



1 - z 



2 1-z) 



1 + z 2 
1 - z 



In z 



+3 + 2z - (J + V) 5(1 - z) 
One therefore arrives at the following formula for the parton cross section 



(46) 



a 9q (z) 



<r (0)— Q{l-z) 

7T 



1 , x / M 2 (l-z) 2 



Q\z 



1 — z 1 

"Cf In z C aa (z) 

2z 2 gqK ' 



(47) 



The g(fci) + g(k 2 ) and (/(fcx) + q(k 2 ) initiated reactions can evidently be deduced 
from the same partonic cross section, as well as the corresponding antiquark 
processes. All these will be included in our numerical analysis. 

2.4 Quark-antiquark annihilation 

Quark-antiquark annihilation (fig. 2a) can be obtained in a straightforward man- 
ner by crossing from the previous calculation: 



\M qq \ 2 (s, t, u) = \M 9q \' 2 (u, t, s) => 32 



991 2/ 



U 



U + tfs 



(46 



No infrared divergency arises and hence e can be put at zero from the beginning. 
The angular integration is trivial and one finds 



a qq (z) 



a s 32 
CT (0) — — z(l - z) 

7T 11 



(49) 



This completes the evaluation of parton cross sections. 



3 Bound state properties 



The decay of r) t are completely dominated by the single quark decay mode [15, |T7 |. 
Decays into two photons constitute a tiny fraction of all events and the branching 
ratio is to a very good approximation given by 

BR( 71 ) = In- (50) 

1 SQD 
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Up to small boundstate corrections T$qd is given by twice the decay rate of a top 
quark T t . For a complete NLO calculation of T] t production and decay into 77 the 
QCD corrections to the decay rates must be included. In Born approximation 



T t (Born) 
Gprrit 2p 



8V2n m t 



1_ — 



m t 



l+(^ 

m t 



mw 
m t 



- 2 



mwY 
m t ) 



(51) 



QCD corrections lead to a reduction by about 7%: 
The complete formula for / is given in DTBf. For the decay rate into a real W, 



2 n 

r t (Born)(l--^/) 

3 7T 



neglecting its width as well as the mass of the bottom quark, one obtains [16[ 



/ 
^1 



2(l-y) 2 (l + 2y) 



Fn 



7r 2 + 2Lz 2 (y)-2L2 2 (l-y) 



+43/(1 - y - 2y 2 ) In y + 2(1 - y) 2 (5 + Ay) ln(l - y) 
-(l-y)(5 + 9y-6y 2 ) 



(52) 



where y 



m 



w 



/m 2 , which is a good approximation to the full answer. The branch- 



ing ratio into 77 is therefore given to high accuracy by 

1 o|i?(0)| 2 ~ 



Br(Q- 



77) 



r 



-12Q^a 2 



SQD 



M 2 



' vr I 3 3 . 



(53) 



The QCD corrections are negative in numerator and denominator and cancel to 
a large extent in the branching ratio. This completes our list of formulae for r] t 
production and decay. 

As stated in the introduction, the most important ingredients for the numeri- 
cal evaluation are the wave function at the origin and the gluon distribution. The 
former are calculated^ for the two- loop potential Vj • The results for the di- 
mensionless quantity i?(0) 2 /M 3 are displayed in fig. 3a for the potentials with the 

0.109,0.1165 



values of 200, 300 and 500 MeV, corresponding to aj^g(M z ] 
and 0.127 respectively (and corresponding to = 0.132, 0.207 and 0.366 MeV). 



This choice for A|^ covers well the present range a-j^g 



MS 



(M z ) = 0.117 + 0.07 de- 
rived from a global fit to the present data 0. In the subsequent discussions the 
results for Vj with the values of 200 and 500 MeV for A will be considered as 
indicative for the uncertainty from the wave function. The branching ratio of r] t 
into 77 is shown in fig. 3b. 

4 The numerical evaluation of R(0) has been performed with the program BOUNDL [Ell. 
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In principle decays into massive gauge and Higgs bosons could be of interest 
and have been calculated in |L5|, |20], |21|Q- The decay rate is dominated entirely by 
the single quark decay. The branching ratio of f] c into 7Z (and correspondingly for 
HZ, ZZ and WW is therefore given by the branching ratio into 77 as calculated 
before, multiplied by the ratio T 7 ^/r 77 . As shown in fig. 3c these latter ratios 
vary between 0.2 (for 7Z) and 40 (for HZ) in the mass range of interest. 



4 Numerical results 

With these ingredients it is straightforward to arrive at numerical predictions 
for the production cross section. If not stated otherwise, these will be based on 



the parton distribution functions of MT (Morfin-Tung) set Bl [ 22]| . a s is chosen 
accordingly to the A value consistent with the parton distribution functions with 
five flavours. All curves correspond to the DIS factorization scheme. 

To investigate the stability of the predictions in a first step the dependence of 
the results on the renormalization scale fi and the factorization scale Qf will be 
explored. The resonance mass will be kept fixed at 240 GeV and /i and Qf will be 
chosen identically. The result of the variation for the gg induced process alone is 
compared to the LO predictions in fig. 4a, the sum of all contributions in fig. 4b 
for V / 5'=16 and 40 TeV. Varying fi = Qf in the range M/2 to 2M induces an 
uncertainty of about ±8% for SSC and LHC These curves have been obtained 
with Vj(200MeV). In the following discussion, both \i and Qf will be fixed to M. 

The magnitude of the QCD corrected cross section including all subprocesses 
is compared to the Born cross section in fig. 5(a) for y/S= 16 and 40 TeV and (b) 
for VS = 1.8 TeV The size of the corrections amounts to 40-50% for LHC and 
SSC energies, whereas the corrections for Tevatron energies amounts to 80-90%. 

To give an idea about the importance of different partonic processes, their 
relative contribution to the production rate is displayed in fig. 6. As anticipated 
in H rjt production is dominated by gluon fusion with the remaining contributions 
amounting about 10% at [i = Qf = M. This holds true for all quark masses and 
energies of interest. 

To study the dependence of the cross section on the choice of parton distribu- 
tion functions, the predictions based on a variety of parametrizations (MT set SN 



22j, DFLM |2|, MRS set B200 J24[, KMRS set B @, GRV j26|) are compared 
in fig. 7 to those based on MT set Bl. ot s is chosen accordingly to the A value 
consistent with the parton distribution functions with five flavours. The resulting 
difference amounts to less than 25% and is in any case small compared to the 
ambiguity from different choices of A in the potential. 



The result for Tzz in [ 21 1 (eq. 9.1c) should be multiplied by a factor (v^ + a 2 a )/lQ. 



12 



The final predictions for the cross section are shown in fig. 8a for LHC and 
SSC and in fig. 9 for 1.8 TeV. The solid line corresponds to Vj(A^ = 300MeV), 

(5) 

parton distribution functions of MT set Bl and a s derived from Aj-^. The dashed 
and dotted curves enclose the range of predictions resulting from the dominant 
uncertainty in the potential and correspond to Vj(Aj^ = 500MeV) to Vj(A^- = 
200MeV). The corresponding predictions for the cross section times the branching 
ratio into 77 are displayed in fig. 8b. 

The results are significantly below those of || even for an optimistic choice of 
the strong coupling constant. 
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Figure captions 

Fig. 1 (a) Lowest order diagrams for eta t production, 

(b) Diagrams contributing from virtual corrections, 

(c,d) Diagrams contributing to corrections from gluon emission, 

(e,f) Diagrams for ggg — > QQ. 

Fig. 2 Diagrams contributing to corrections from (a) qq — > r) t q and 
(b) qg — > r] t g gluon emission. 

Fig. 3. (a) Predictions for i?(0) 2 /M 3 of the S-wave ground state and 

(b) branching ratio of i] t — > 77 as functions of M for the potential Vj with 
A^L=200 (solid), 300 (dotted) and 500 MeV (dashed). 

(c) Decay rates of rj t into 7Z, ZZ WW and HZ (for m H = 70 GeV) 
normalized to the 77 channel. 
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Fig 4. (a) Scale dependence for pp —>■ rj t + X for the gg induced process alone 
and (b) scale dependence for all contributions. Both LO (dashed) and 
NLO (solid) results are shown for v^=16 and 40 TeV. We use MT set Bl 

DIS factorization scheme. 



parton distributions with A-=-r ?T =194 MeV for four flavours and work in the 

1 MS 



Fig 5. (a) Ratio between the radiatively corrected cross section for pp — > i] t + X 
and the lowest order result for \/S=16 and 40 TeV and (b) ratio between 
the radiatively corrected cross section for pp — > rjt + X and the lowest 
order result for V~S=1.8 TeV. We use MT set Bl parton distributions with 
A^^=194 MeV for four flavours and work in the DIS factorization scheme. 

Fig 6. Relative contributions to rjt production at /i = Qp = M for "vS—16 

(dotted) and 40 (solid) TeV. LO contributions are (A) gg —*■ r] t + X, NLO 
contributions are from (B) gg — > r\ t + g + X , (C) gq — > r\ t + q + X , (D) 
(qq -> 7ft + 5 + X) x 100. 

Fig 7. Prediction for the production cross section of r] t at LHC using different 
sets of parton distributions functions (MT set SN @, DFLM @, MRS 



set B200 [g|], KMRS set B |2§ , GRV@) normalized to the predictions 
obtained with MT set Bl " V1 



Fig 8. (a) Cross section for rjt production including NLO corrections and (b) 
cross section multiplied by the two photon branching ratio at \/~S = 16 
and 40 TeV for the potential Vj with A^L=300 (solid), A^=200 (dotted) 
and A^=500 MeV (dashed). 

Fig 9. Same as fig. 8a for 1.8 TeV. 
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